On the temporal Wilson loop in the Hamiltonian approach in Coulomb gauge 



We investigate the temporal Wilson loop using the Hamiltonian approach to Yang-Mills theory. 
In simple cases such as the Abelian theory or the non-Abelian theory in (1 + 1) dimensions, the 
known results can be reproduced using unitary transformations to take care of time evolution. We 
show how Coulomb gauge can be used for an alternative solution if the exact ground state wave 
functional is known. In the most interesting case of Yang-Mills theory in (3 -I- 1) dimensions, the 
vacuum wave functional is not known, but recent variational approaches in Coulomb gauge give 
a decent approximation. We use this formulation to compute the temporal Wilson loop and find 
that the Wilson and Coulomb string tension agree within our approximation scheme. Possible 
improvements of these findings are briefly discussed. 
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I. INTRODUCTION 



In recent years there has been a renewed interest in studying Yang-Mihs theory in Coulomb gauge, both 
in the continuum [ll and on the lattice These developments are mainly based on the fact that 

the consequences of the Gribov-Zwanziger picture of confinement 0, [lo| are particularly transparent in 
Coulomb gauge: Physical states in the Hamiltonian formulation have to obey Gauss' law, without the 
need of additional restrictions such as a vanishing color charge based, in turn, on the assumption of a 
globally conserved BRST charge. Furthermore, since Gauss' law can be resolved explicitly in Coulomb 
gauge, any renormalizable Ansatz for a vacuum wave functional is admissible in this gauge without 
having to explicitly construct the physical Hilbert space [ll|, which is a key fact to make variational 
approaches viable. (The equivalent Gupta-Bleuler type of constraints in covariant gauges have no such 
simple resolution.) 

Consequently, much work was devoted to a variational solution of the Yang-Mills Schrodinger equation 



in Coulomb gauge [12H20|. Using Gaussian type of Ansatze for the vacuum wave functional, a set of 
Dyson- Schwinger equations for the gluon and ghost propagators was derived by minimizing the vacuum 
energy density. An infrared analysis of these equations exhibits solutions in accord with the Gribov- 
Zwanziger confinement scenario. Imposing Zwanziger's horizon condition one finds an infrared diverging 
gluon energy and a linear rising static quark (Coulomb) potential [l^ and also a perimeter law for 
the 't Hooft loop [20]. These infrared properties are reproduced by a full numerical solution of the 
Dyson- Schwinger equations over the entire momentum regime (l6| . and are also supported by lattice 



c.,c.,.t.o„. U^. M„.„ve,, Since the ™ of the gho.t d,e.i„g fa„ct.„„ can be idenUfled with the 

dielectric function of the Yang-Mills vacuum 2l|, the latter behaves in the infra-red, by the horizon 
condition, like a perfect color dia-electric medium, i.e. a dual superconductor. This suggests that the 
linearly rising Coulomb potential is a manifestation of confinement via the dual Meissner effect. 
These nice physical conclusions are, however, a bit premature. Using variational arguments, one can show 
that the Coulomb potential is only an upper bound for the true potential between heavy quarks, which 
implies "no confinement without Coulomb confinement" but not the converse. To get a sufficient 
criterion for confinement, one has to study the true potential extracted from the temporal Wilson loop. In 
the non-Abelian continuum theory, the Wilson loop is difficult to calculate because of path ordering. In 
ref. (221, the spatial Wilson loop was computed from a Dyson-Schwinger equation, which has previously 



been derived in ref. [23] in the context of supersymmetric Yang-Mills theory. Although a linearly rising 
potential was extracted from the obtained Wilson loop, it is not clear to which extend this approach also 
applies to non-supersymmetric Yang-Mills theory. 

As was clarified in ref. [2J], it is by no means trivial to extract the static potential from the Wilson loop 



even in simple cases such as QED, due to the presence of singular self-energy terms. Another attempt 
to extract the static Coulomb potential of QED from the Wilson loop was undertaken in ref. where 
a partial wave expansion has been used. We have not been able to reproduce the result of ref. [25| for 
the static potential in QED, and we reconsider the relevant partial wave expansion in appendix [A] In 
the present paper, we will instead show that the proper and elegant way to extract the QED potential 
from the Wilson loop should be based on a separation of longitudinal and transversal degrees of freedom. 
In Coulomb gauge, this separation is accomplished automatically^ which is the main reason why this 
particular gauge is so convenient for our purposes. 

In the present paper, we will attempt to extract the static quark potential from the Wilson loop within the 



FIG. 1. Temporal rectangular Wilson loop. 



Hamiltonian approach, using gauge-invariant and Coulomb gauge techniques in exactly solvable cases, and 
approximate results from variational calculations for the most interesting case of non-Abelian Yang-Mills 
theory in D = 3 + 1. 

Our paper is organized as follows: In the next section, we demonstrate how a unitary transformation 
similar to the one that leads to the interaction picture in quantum mechanics gives rise to an induced 
electric field that allows to compute the Wilson potential exactly in special cases. To prepare for the 
non-Abelian theory in (3 + 1) dimensions, we also show how Coulomb gauge can be incorporated in these 
exactly solvable models. In section 3, we first treat the case of an Abelian U{1) theory reproducing the 
usual Coulomb potential from Maxwell's theory. Section 4 treats Yang-Mills theory in (1 + 1) dimensions 
on a cylindrical spacetime manifold. Finally, section five employs the known Gaussian wave functional 
from variational approaches to compute, with additional approximations, the Wilson string tension, which 
happens to agree with the Coulomb string tension within the given approximation scheme. We conclude 
with a brief summary and an outlook on further improvements of our calculation. 

II. THE TEMPORAL WILSON LOOP IN THE HAMILTON APPROACH 

Below, we investigate the temporal Wilson loop in the Hamiltonian approach to SU{N) Yang-Mills 
theory. This formulation assumes the Weyl gauge, Aq = 0, leaving a residual symmetry under spatial 
gauge rotations. Since the Wilson loop operator is gauge invariant, it is not necessary to further fix this 
residual symmetry in all cases. In fact, we will show that a gauge-invariant Hamiltonian treatment of the 
temporal Wilson loop is successful in simple cases, while more complex situations require gauge fixing as 
a technical tool to simplify the separation of physical degrees of freedom. 

A. Physical interpretation 

We consider a planar time-like Wilson loop as shown in figure [T] If the (Euclidean) time coordinate 
r = xo is unrestricted, the Weyl gauge Aq = can be adopted and the Wilson loop operator in the quark 
representation r reads 

?7[A](R,0; T)1" • J7[A](R,0; 0)1 , (1) 



W = —tTr 



[/[A](0,R; T = T) • [/[A](R,0; r = 0) 



tr 
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where ir^ is the trace in the SU{N) representation r, dr = tr,. 1 is its dimension and tr = d^^ tr denotes 
the normahsed trace. Furthermore, the parallel transporter along the path running from xi to X2 at 
the time xq = t is given by 



U[A] (x2,xi ; r) = P exp 



J dxA(x, r) 



(2) 



In this expression, we have absorbed the gauge coupling g in the connection A, which entails that the 
Yang-Mills Hamiltonian (in the absence of gauge fixing) takes the form 

H = ^ [ dxn^ + ^ [ dxB^ , (3) 



2 J 2g 

The time evolution of the gauge field in Euclidean space is then formally 

A(x,r) = e-^^ A(x,0) e"^ . 

and it is easy to see that the same semi-group law applies to the parallel transporter, so that the vev of 
the Wilson loop becomes 

( ) = tr ( I f/t (R) e-(^-^°)^ f/(R) I ) = tr ( R | e'^^^^^") | R ) . (5) 

Here, Eq is the energy of the YM vacuum |0) and the abbreviation C/(R) ~ C/[A](R, 0;r = 0) for the 
straight line transporter — > R made from gauge fields at time t ~ was introduced. 

The next step is to insert a complete set of energy eigenstates in the matrix element eq. ([5]). To do 
this properly, we have to recall that the Hamiltonian eq. ([3]) is accompanied by a Hilbert space which 
decomposes into charge subspaces by virtue of Gauss' law, 

f ''(z) (A ; XI, . . . , x„) = - I E E [^{re}] ^ ^ ^(^ - n \ ■ ; XI, . . . x„) 

(6) 

where the sum is over the color spins of static quarks with color representation Tr located at positions 
xi,...,x„. Since the Gauss law operator r''(x) = D°''(x)n''(x) commutes with the (gauge-invariant) 
Hamiltonian, the charge subspaces are invariant under the YM dynamics and mutually orthogonal. By 
Schur's lemma, the Hamilton operator eq. ^ is thus proportional to the color unit matrix when acting 
on any of the charge subspaces. Moreover, the Gauss law operator generates time-independent gauge 
transformations (7(x), i.e. the iteration of eq. ([6|) states (in the absence of a vacuum angle 6 = 0) that a 
physical state from the n-quark sector transforms under gauge rotations as 

*,„...,,„(A^;xi,...,x„)= E 47-![t^(xi)]---<'fc![t^(x«)]*fci,...,fc.(A;xi,...,x„), (7) 

ki ,....kji 

where f^*^*"^ [fl] denotes the color matrix fl in representation r. For notational simplicity, we will henceforth 
assume that the Wilson loop is in the fundamental representation, so that all representation symbols i^'-''-' 
can be omitted. 



It is now trivial to see that the Wilson state [26|, |27| |R) = ?7(R) |0) from eq. ([5|) corresponds to the 
charge sector with a quark at position R and an antiquark in the origin. In fact, the identity 

r"(z)[/(R) = -S{z - R) iT" U(R) + 6{z - 0) U{R) iT" , (8) 

implies the matrix relation 

r'^(z) I R),, = ~{,5(z - R) {tT%k 5,1 - (5(z) 5,k (*r«),,} | R),., , (9) 

which corresponds to eq. ([6|) . For finite gauge transformations 51 acting on the Wilson state, the equivalent 
of eq. d?]) is 

|R)„- ^ Sl(R),fc • |R)fc£ • 52(0)7/. 

Finally, for quark and antiquark in the same color representation, it is often convenient to parallel 
transport the charges to the same local color frame so that the charge density matrix acts on a single 
spin index only. In the present case, this can be achieved by the identity 

[U-i^T'^)]u-[U^T'^U\^-UM (10) 

so that eq. ^ takes the simpler form 

r'^(z) {/,,(R) = {(5(z) [;7(R) iT- [/(R)t],fe - 5(z - R) zT^} {/^.(R) = -p^,(z) • Uu{B.) , (11) 

where the second spin index [l] is not affected. As a consequence of the parallel transport, the two 
charges at the endpoints of the Wilson line do not exactly compensate 

= y"dz/9''(z) = i [C/(R)r''C/(R)t -T'^] ^ (12) 

unless the color group is Abelian. (In any case, however, trQ" = 0.) 

Returning to eq. (O, we insert a complete set of energy eigenstates from all charge sectors, but only the 
states from the proper quark-antiquark sector will give a non-vanishing contribution, 

n n 

where = tr(f2'l'57). In the limit of large Euclidean time extensions, we project out the static quark 

potential 

- lim T-^ In IW) = i;^'^'^^ -Eq = V(R) , 

provided that the Wilson line has non-vanishing overlap with the true vacuum in the proper quark- 
antiquark sector. (We will come back to this issue in section HVl below.) 

B. Induced electric field 

To proceed, we use the spatial parallel transporter U (R) to induce a transformation similar to the one 
that leads to the interaction picture in quantum mechanics. For instance, the new Hamilton operator 
becomes 

H = U\R)HU{Il) = H(R). (14) 
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All transformed quantities will be R-dependent and denoted, generally, by a tilde; they will also be color 
matrices in the representation chosen for the Wilson loop, although we will sometimes omit color indices 
for simplicity. For the transformed momentum operator, we obtain 

n^(x) = ;7t(R) n^(x) c/(r) = n^(x) + e^(x) , (15) 

where a new color electric flux field 

e'^(x)^ J dy6{^-y)zU{yyT^U{y)= J dy <5(x - y) "^(y) (16) 

with T°'{y) = [/(y)^ [/(y) has emerged. The unitary transformation (IT5t removes the parallel trans- 
porters from the Wilson state so that the Wilson loop becomes the expectation value of the transformed 
time evolution operator in the Yang-Mills (zero charge) vacuum, 

(M^) =tr(0|e-^(^-^°) |0). (17) 

The effect of the Wilson lines is now shuffled into the kinetic term of the new Hamiltonian, 



j [n''(x) 1 + e"(x)]' + ^ j Sx-B{^f 1 . (18) 

Note that e''(x) depends on the gauge field, so that [11° (x), e''(x)] ^ which complicates the treatment of 
the new Hamiltonian considerably. The induced electric field satisfies Gauss' law for the external charge 
density induced by the Wilson line, 

[D'^''(x), e''(x)] = J7(R) V(x) C/(R) = iT"^ (5(x) - ^f '^(R) (5(x - R) = p'^(x) . (19) 

This charge density p(x) is again a color matrix similar to eq. (jlip . 

Before proceeding, let us briefly discuss the piece in the transformed Hamiltonian (llSp which is quadratic 
in the induced electric field, 

V,^Am = ^^ j d'x[e'^{^)Y . (20) 

After a simple calculation using the definition of the quadratic Casimir operator T"^ = T° T° = C2 1, 
we find 



52 



J J dyS\^-y)l. (21) 



This term is independent of the gauge field and thus also of the vacuum wave functional used, but it does 
depend on the shape of the contour : R in the spatial part of the Wilson loop. Choosing to 
be a straight line connecting and R, we have 

^ndi"^] - yC2<5i(0)-|R|l . (22) 



The color unit matrix disappears when the trace of the Wilson loop is taken. While eq. ([22| is formally 
a linearly rising potential, the corresponding string tension 

3lc,6liO)^^g^C2 (23) 
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is UV-divergent (cutoff A) and of pureiy fcinematic origin. Except for tlie Casimir operator C2 — > 1, the 
same contribution is also found in the case of QED which does not confine; so it must be spurious. In 
fact, the singular energy is due to the infinitely thin Wilson lines whose point-like cross section gives 
rise to the factor 5^(0) in eq. (|23p . If the lines are smeared out to an effective tube with cross section 
^ A^^, the self-energy is regularized, (5^(0) ^ A^, but it will still obscure the physical potential unless it 
is properly isolated and removed. 

To see the origin of the problem in more detail, consider the Abelian case where Gauss' law for the induced 
electric field, Ve ~ p, indicates that only the longitudinal part of e is sensitive to external charges. As a 
consequence, the static potential (|20l) decomposes into two parts. 




The first piece V± is built from alone and thus insensitive to external charges; it represents a vacuum 
property that renormalizes the Wilson loop but does not affect the physical potential between static 
quarks. The latter is given exclusively by the second piece V\\ which will be shown below to describe the 
well-known Coulomb interaction. The separation of transversal and longitudinal degrees of freedom is 
thus essential to disentangle the physical potential from renormalization effects. In rcf. [25^, a different 
approach to remove the self-energies of the Wilson line and to isolate the proper quark potential was 
proposed, which is based on partial wave expansion. Appendix lAl compares this method to our reasoning 
based on Gauss' law. 



C. Coulomb gauge 

Gauge invariant calculations of the Wilson loop are only possible in special cases where the exact ground 
state is known analytically. In realistic cases, the vacuum wave functional is unknown and one has to 
resort to approximate calculations which usually require some sort of gauge fixing. For the Hamiltonian 
approach. Coulomb gauge is a natural choice which simplifies the separation of dynamical from gauge 
degrees of freedom. In QED, for instance, the Coulomb gauge fixed field represents the physical (gauge 
invariant) degrees of freedom. 

Not surprisingly, the variational approach to YM theory in Coulomb gauge has recently become an 
excellent tool to obtain approximate ground state wave functionals which are suitable for practical com- 
putations. We will use this approach for YM-theory in 3 -I- 1 in section IVl 

In view of our previous analysis, we can choose to fix the Coulomb gauge either before or after the 
unitary transformation that led from expression ([5]) to eq. p?| . That is, starting from the invariant 
Hamiltonian H in eq. we can either 

1. gauge fix H directly and compute (W) from eq. ([5]): (W) = tr ( R | e-T{H-Ea) | j^^ 

The gauge fixed Coulomb Hamiltonian can be obtained either by canonical or by path integral 
techniques. The well-known result is the expression [^l 

^=y 1^3x^-1 n^(x)^n^(x) + ^ j d'x[B%^)]^ + Hc (24) 
^C = Y J ^'(^'2^) ^°(^) + '°Sy„(x)] ^ • F'"'(x,y) [p\y) + p^y„(y)] , (25) 
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where ^[A\ = dct(— DV) is the Faddeev-Popov determinant and the Coulomb kernel F reads 

F = (-DV)-i(-A) (-DV)~i . 

The resolution of Gauss' law must take into account that the state |R) in eq. ([5|) is charged, 
i.e. besides the dynamical charge of the gluons, /o3y„(x) = -Af" (x) n'^(x) (with Af" = f""^ A'^^), 
we also have to include the external charge induced by the Wilson line. For instance, from Gauss' 
law in the form eq. (fTT|) 

the usual resolution entails that the charges in eq. are really matrices [p° + pj^^^^ 1] . and so is 
the Coulomb kernel He- The Wilson state | R), after gauge fixing, depends on A-'- only and is 
thus no longer from the gg-sector; it will, in fact, have overlap with the zero-charge sector. 

2. compute H = f7(R) H C/(R)t and gauge fix expression eq. (H?]): (ly) = ( | tr e-'^^^-^") \ ) 

Since the Wilson lines have been taken care of by the unitary transformation iJ — > in eq. , 
no further external charges have to be included, and Gauss' law takes the form DII | ) =0. All 
calculations can thus be done in the zero-charge sector. From eq. (I18p . we could now follow the 
standard route and apply the Faddeev-Popov procedure to the vev 

(0|il|0) =y"^A[(n + e)*o]* [(n + e)*o] +••• 1'o[A] = (A|0). 

Using Gauss' law to solve for n" eliminates all longitudinal d.o.f. and yields the vev ( | \ ) 
in terms of transversal fields only. However, this is usually not sufficient to find the required vev of 
the gauge fixed evolution operator (0 \e^'^^ | 0), except for special cases such as QED. 

Thus for the non-Abelian case, we prefer to follow the first method based on the standard Coulomb gauge 
Hamiltonian, while Abelian systems can be treated with both approaches. 



III. QUANTUM ELECTRODYNAMICS 



A. Gauge invariant treatment 



In the Abelian case G — U{1)^ the Maxwell gauge potentials do not carry a color index, and the hermitian 
generator of the non-Abelian gauge group, (jT°), has to be replaced by 1, so that also C2 = 1. The formula 
for the induced electric field in the representation eq. ([18]) simplifies to 



(26) 



e(x) = j dy5(3)(x-y). 
In particular, e(x) is field-independent and the parallel transporter — > R becomes 



C/(R) = exp 



i J d^x e(x) • A(x) 



To isolate the physical degrees of freedom, we split all vector fields in longitudinal and transversal parts. 
Without gauge fixing, the ground state vJ/q [A] — ( A | ) is eigenstate of the gauge invariant Hamiltonian, 
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and thus gauge invariant itself; as a consequence, it can only depend on the physical (i.e. transversal) 
field , so that n" | 0) = 0. From eq. ^ and this entails 



(W^) = 



-^QED ~ 2 



(27) 



Since e is independent of the gauge field, the last term in ^Jqeq commutes with the remainder and can 
be pulled out of the vev, 



'W) =e- 



-TV" 



g-T(ffQED--Eo) 

1^11 = ^ |d3x[ell(x)]^ 



(28) 



where Vn is the part of the induced static potential that is sensitive to external charges, cf. our previous 
discussion. The remaining Hamiltonian 



1 



[n^(x) + e^(x)]^ + -[VxA^(x)] 



9^ 



depends on the physical degrees of freedom only. It is most easily treated by reversing the transfor- 
mation eq. p^. but this time with the parallel transporter made of physical gauge fields A^ only, 



C/^(R) EE exp 



I I dxA (x) 



= exp 



ij d3xe-L(x) • A^(x) 



(29) 



As a result, the induced transversal electric field is removed from the Hamiltonian and re-shuffled into 
the state | ) = [/-^(R) | ). Although this new state | ) resembles the Wilson state, it depends on 
A^ only and is thus gauge invariant, i.e. | R^ ) lies in the zero-charge sector. Inserting a complete set of 
zero-charge (gauge invariant) eigenstates of the standard QED Hamiltonian, 



^QED I d"x 



1 



[n^(x)]^ + -[Vx A^(x)] 



with = 5/i6A-'- we find 



[W)=e-^^ii ^e-(^"-^«)^ |(n |C/^(R)| 0) 



(30) 



(31) 



It is clear that only the first term (n — 0) in the sum survives in the limit T — )• oo, provided that the 
matrix element ( | ;7-L(R) | ) = ( | R^ ) is non- vanishing. In the present case, this matrix element can 
be computed explicitly, since the vacuum wave functional is Gaussian, 



A-L(x) 

A simple Gaussian integral gives 







^ .vK exp 







dHx,y) A^(x)c.(x,y) A^(y) 



;0|C/-(R)|0)= 



exp 



; J d^xA^ix) ■ e-L(x) 
I d3(x,y)e^(x)c.-i(x,y)e^(y) 



= 



exp 



R 

/ dx,. 



R 



dy-jUj{x.)uj ^(x,y) 
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where tij is the usual transversal projector. This expression can be further evaluated if we choose the path 
^ R to be a straight line and employ the QED dispersion relation w(k) = |k| to invert the kernel w(x, y) 
in momentum space. Introducing spherical coordinates in momentum space, the parameter integrals for 
the Wilson lines can be done first, while the subsequent angle integral can be expressed through the 
integral sinus, Si(x). For the remaining integral over k = |k|, we use a sharp UV cutoff (i.e. a 0(3) 
invariant cutoff in 3-space), to find 

( I ) = exp ("47 + 2 cos(Aii') - 4 In(Aii') + 4 Ci(Ai?) - 2 



AR 



2Ai?Si(Ai?)j 



exp 



In the limit T 



where 



UAR - 4 ln(Ai?) - 47 + ^(A-3)) 
00, the bare Wilson loop thus becomes 

{W} = Zw (AR) ■ exp [-T Vj, (i?)] 



(32) 



(33) 



Zw{x) 



exp 



9 ( 



4 In X - 47) 



Formally, Zw vanishes in the limit A 00, which indicates that the state | R^ ) with support on 
an infinitely thin line has poor overlap with the true QED ground state. On the other hand, Zw 
is independent of T and therefore does not contribute to the potential. The correct interpretation is 
given by the operator product expansion (OPE) [29]: Since the Wilson loop ly is a composite operator 
involving products of arbitrarily many gauge fields, we expect short-distance OPE divergences associated 
with the operator W ^ in addition to those counter terms which are necessary to render elementary Green 
functions finite. Eq. (|33|) shows that this is indeed the case, with one overall multiplication of the operator 
W — >■ Wu = ZZ} ■ W rendering all matrix elements finite. This renormalization property of Wilson loops 



is an exact result that has been known for a long time both in the Abelian and non-Abelian case 
From the renormalized value 



30 



31|. 



; Wfl ) = exp 



-TV«{R) 



we conclude that, in the Abelian case, the Wilson potential agrees with Vj| . A simple explicit calculation 
gives first the induced electric field 



e"(x) = -V,(-A,i)V,e(x) = V,(-A,i) / dyVySix-y 

Jo 



-V, (~A-i) h<5(x - R) + <5(x)] ^ -V, (-A-1) p(x) 



(34) 



and thus Gauss' law Ve" — p. The Wilson potential becomes, after integration by parts and using the 
Green function G(x, y) = ( x | (- A)"i | y ) = (47r|x - y|)^^ 



ell(x) 



' = ^ / d3(x,y,z) -V,G(x,y)p(y) • - V, G(x, z) p(z) 



I d3(x,y)p(x)G(x,y)p(y) = ^ J d^{^,y) 



p(x) pjy) 

Att |x - y| 
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9 



47r|R| 47r|0| 

This is just the familiar Coulomb law including the self-energy of the two charges. 



(35) 
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B. Coulomb gauge 

In the Abelian case, we can apply both methods described in section Hi CI to evaluate the Wilson loop in 
Coulomb gauge: 

1. In the Abelian case, the gauge fixed Hamiltonian eq. ([M)) simplifies considerably. From ^ = const, 
F = (— A)~^ and pdyn = 0, we have 

^2 



^^2 



i jd?x g2n^(x)2 + lB(x)2 +|- Jd^x,y)p{^)i-A)-'i^,y)p(y)^HQEB + Hc. 



The second piece He is a field- independent c- number, so that e~^^ = ^-tHqed . q-tHc _ Taking 
into account that | R) = ?7[A-'-] | ) and inserting a complete set of eigenstates of Hqed, we arrive 
directly at eq. (pij) . from which we can follow all the steps to the final expression eq. ([55)) for the 
Wilson loop. 

2. Resolution of Gauss' law in the chargeless sector gives II" | 0) =0, and after splitting into longitu- 
dinal and transversal fields and Coulomb gauge fixing, the operator H in eq. (|18p becomes i?QED 
eq. (P7| . From there, we can again follow all the steps that led eventually to the expression eq. (pSj) 
for the Wilson loop. 



IV. YANG-MILLS THEORY IN D = 1 + 1 

In the non- Abelian case, the exact vacuum wave functional is only known for the special case of one space 
dimension, which is nonetheless a good testing ground since the exact result for the string tension is known 



e.g. from path integral techniques [32|, |33|. However, YM theory in _D = 1 + 1 is (almost) topological 
and non-trivial results can only be obtained on non-contractible Euclidean spacetime-manifolds. Our 
Hamiltonian approach at zero temperature requires that Euclidean time be unrestricted, so that the only 
non-trivial spacetime manifold is M = R x S*^, where the first factor is time and the second factor is 
space. In the following, we will therefore assume that the spacetime is a cylinder, i.e. space is the interval 
[0, L] with periodic boundary conditions. (The general theory of fiber bundles indicates that transition 
functions on M can be taken trivial, i.e. gauge fields are strictly periodic while gauge transformations are 
periodic up to a center element.) 



A. Gauge invariant treatment 

In D = 1 -I- 1, there is no magnetic field and the vacuum is characterized by the condition 11° (a;) | ) = 0. 
This is because then {Q \ H \ 0) = which in view of 77 > is an absolute minimum so that | ) must be 
the ground state with zero energy Eq = 0. The condition H*^ | ) = can be readily solved by a constant 

*o[^] = (^|0) const. (36) 

If we first ignore the commutator [11° (x) , e'^ (x)] 7^ in eq. , we find immediately H \0) = Vi^d [^] '10), 
where Vindi*^] is the static potential defined in eq. pi[) . As a consequence, we have 

(M^) =tr(0|e-^^|0) =tre-^^'-['^l 
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FIG. 2. Two possible Wilson lines connecting the same points on the compact space interval with periodic 
boundary conditions. 



i.e. the Wilson potential matches the induced potential eq. (|2ip . 

"1 i-R i-R. 



V{R) = Vn,^[0^R] = ^C2l f dx I dy6{x,y)^^C2Rl. (37) 
^ Ja Jo ^ 



10 Jo 

We are thus led to a linearly rising potential with the string tension 

-,2 



^1+1 = Y . (38) 



This is the correct result known e.g. from path integral techniques 32ll33l| or from the lattice Hamiltonian 



ll|. However, there are two caveats to this result: 

• we have neglected the commutator [H'^ {x) , e'^ (x)] in the Hamiltonian ; 

• we have not yet made use of the periodic boundary conditions and the resulting potential eq. p7p 
is not invariant under the replacement R ^ L — R, nor is it periodic in R. 

We will not try to solve the first problem here, since the correct treatment of the dynamics will be 
easier in Coulomb gauge as described in the next section. As for the periodicity, the distances R and 
L — R are indeed equivalent on a spacetime cylinder, but the Wilson loop will not be invariant under 
the replacement i? — >■ L — R. This is not a contradiction since the physical potential in the compact 
spacetime setting is no longer given by the exponential of any single Wilson loop. 

To understand this subtle point in more detail, we have to go back to the physical interpretation of 
the Wilson loop in section IH Al The identification of the (exponent of the) Wilson loop with the static 
quark potential relies on the assumption that the Wilson state made by the direct path : — > i? has 
non-vanishing overlap with the true YM ground state in the qq-sectov for all values of R. However, the 
compact space interval allows for many more Wilson states from the same gg-subspace, namely (i) the 
state where the path 'rf winds around the compact space in the opposite direction (cf. figure [2]), and 
(ii) all states which can be obtained by winding the Wilson line an additional integer number of times 
around the space in any direction. We can label all these states by the end point of the Wilson line 
used to create them, || i? + mL) \m ^TL^. (A negative end point indicates that the Wilson line winds 
around space in negative direction; for instance, the two lines in figure [2] correspond to the states | R) 
and I i? — L ), respectively.) 

At small i? <C we expect that the state is energetically favored, but this may change at larger 
R so that the Wilson loop with the direct path connecting R will no longer give the correct quark 
potential. Assuming that at least one of the Wilson states | R + mL ) has overlap with the true qq ground 
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state, we can extract the true quark potential from the largest eigenvalue of the matrix 

(i? + mL|e-'^(^~-^o) |i? + nL), m, n e Z . (39) 

The same information may also be obtained from the diagonal matrix elements, 

ViR) = min [ - lini T^^ \n{R + mL\ g-^t^-^o) \R + mL)] . (40) 

This follows from section Hi Al because the quantity in square brackets is simply {E^^^f^-^ — £^o)i where 
^e(m) I'fifers to the lowest Yang-Mills eigenstate in the gg-sector which has non-vanishing overlap with 
the Wilson state | R + mL ) . Assuming that at least one of these Wilson states has overlap with the 
true qq-vacuum, we must have €(rn) = for at least one m. Thus taking the minimum with respect 
to m yields (Elf^^ - Eq) on the rhs of eq. (PH)) . which is the proper static quark potential. However, 
since the minimum can (and will be) be at different m for different distances R, the true quark potential 
is not obtained from any single Wilson loop. It is already evident that the potential (PU)) has all the 
required periodicity properties; the explicit evaluation of the matrix elements is, however, best performed 
in Coulomb gauge. 



B. Coulomb gauge 

Our treatment of Yang-Mills theory in _D — 1 + 1 in Coulomb gauge will follow the first method outlined 
in section HlCi i.e. we gauge fix the usual Yang-Mills Hamiltonian directly, before removing the Wilson 
line from the states. In (1 -I- 1) dimensions, the Coulomb condition 

diAf = (41) 

leaves a (spatially) constant gauge field — const, which can be further diagonalized by means of the 
residual global color symmetry [diagonal Coulomb gauge) [isl. [s^. For simplicity, we will restrict the 
color group in this subsection to G = SU{2) with the antihermitean generators T° = T°/(2i), where 
are the Pauli matrices. In the diagonal Coulomb gauge, we are then left with only one physical degree of 
freedom, A^^^, which we normalize to the dimensionless variable 

i^ = ^AlL€[0,n]. (42) 

The restriction to the compact range — the fundamental modular region — ensures uniqueness of the 
gf. section and avoids the Gribov problem. Using this variable, the gauge-fixed Hamiltonian in the absence 
of external charges becomes (isl 

i f' qL\^ 1 d ^ „ d , ^ 

H = — — ^ s n^ . 43) 

2L\2 ) sin^ ^ d-d d-d ' 

The corresponding eigenfunctions H\n) — En\n) are just the characters of SU{2) and thus known 
explicitly, 

, , , [2 sin(n + l)i? 1 f qLV , , 
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In diagonal Coulomb gauge, both A and II live in the Cartan algebra, so that their color commutator (the 
gluon charge) pdyn vanishes. As a consequence, the Coulomb term He from eq. ((25|) vanishes completely 
unless external charges p ^ are present. For the Wilson loop, the external charge p is given by eq. pT|) : 
it is field dependent but contains no derivative operator II and thus commutes with the Jacobian ^ . 
We are thus left with 

hc^yJ p^^""^ '"'(^^ • ^^^^ 

Notice that this expression differs from the static Coulomb potential because the charges p{x) from 

eq. (jlip still contain the parallel transporters from the original Wilson lines. 

To further compute eq. (1451) . it is convenient to introduce a polar basis in color space, 

/I 



ecr=l 



1 

71 



















4 - 


(e-)*r„ = (-l)- 



(46) 



and expand all relevant color vectors in this basis. 

Pa = el Pa 

(e^)*F-*(x,2/)e^ = b,,F„{x-y). 
After a simple calculation, we find the rotated generators 

where \J is the parallel transporter ^ i? + mL in diagonal Coulomb gauge. 



(47) 



U — exp (ia T3) 



a 



R 1 
— + m I -d . 



(48) 



We can now insert this representation in eq. (|lll) to find the external charge of the Wilson state | R + mL ] 
in the form 



o,{x) = -S{x -R-mL) + 5{x) e2'-(f+™)'' , 



(49) 



As mentioned earlier, this differs from the sum of two static point charges by the Wilson line connecting 
them. The Coulomb Hamiltonian eq. (|45p can now be put in the form 

.2 +1 



Hr 



9 



rt 



Hciz) 



cr=-l 



(50) 



where 2 = i? + mL. To complete the calculation of He, we can use the explicit expressions for the 

( — I 

Coulomb kernel derived in ref. [18] to find the polar basis representation 



L 

2^ 



E 



cos(27rnz/L) \z\ 
^ 2Z ^ Y 



L 
12 



r 00 



^±(^) = z E 



{■KTl ± §Y 



L \z\ z 
T TT T^t: coi^ 



(51) 



0.25 0.5 0.75 I 0.25 0.5 0.75 1 

R/L R/L 

FIG. 3. Left panel: The Coulomb Hamiltonian Vc{z) for the four shortest Wilson states 1^). Right panel: The 
true quark-antiquark potential for D = 1 + 1 Yang-Mills theory on the spacetime cylinder, as extracted from the 
prescription eq. (|40[) . 



The Fourier representation shows that the Coulomb kernel is indeed periodic, F„{R + mL) = Fa^{R), 
while the explicit expressions on the rhs. are only valid for \z\ < R, with periodic continuation outside 
this range. After a brief calculation, we find, again in the range \z\ < L, 



hoiz) = \z\ 



L ' 



h±{z) = \z\ . 



With To-rJ = 1 + (TTa, the Coulomb term is finally given by 



■,2 r 
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3UI 



\z\<L. 



(52) 



(53) 



Notice that this formula is only valid for \z\ = |i? + mL\ < L, i.e. for the two states | R) and \ R ~ L) 
from figure [2j For other Wilson states, we have to go back to the Fourier representation in eq. ((5T|) 
and do the sums numerically. It is then easy to see that for < i? < L, only the two states | R) and 
\R — L) compete for the minimum, while all other Wilson states have higher energy, cf. figure [31 Since 
the Coulomb Hamiltonian is color diagonal and field-independent for these states, we can complete the 
calculation of the relevant matrix elements as follows: 



tr(z|e-^(^+^^)|z) =tre-^^^(^)• 



^|e-^^U; 



OO 

n=0 



-^°)tr \{n\U{z)\0)\^ , (54) 



where z = Rorz = R — L and Hc{z) — Vc{z) 1. In the last step, we have inserted a complete set of 
eigenstates from the zero quark sector, since the Wilson state \z) made from transversal gauge fields is 
no longer in the qq sector and expected to have overlap with the zero-charge ground state. In fact, the 
overlap can be worked out explicitly with the result 



tr 



(n|[/(z)|0)|' = -|(l)'(l_(„)"cos|^) 



1 



2f 



(55) 



The final sum over energy eigenvalues n in eq. H54p can only be performed numerically. At large spacetime 
volumes LT ^ .9^^, the ground state n ~ dominates the sum, and all remaining contributions are 
exponentially suppressed. The relevant matrix element |( | i? ) p is then independent of T and approaches 
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unity for large space extensions L ^ R. This entails that 

z e{R,R-L}. 
1 competing for the minimum ii \R\ < L) yields the 



- lim T-'^ \ntT{z\e~'^'^"-^"^z) = Vc{z) 

T— !-oo 



The prescription (|40p (with only m — and m 
final result 



V(R) — min VciR + rnL)— min 

me{o,-i} me{o,-i} 



3 o , , / IR + mL 



(56) 



This function is plotted in the right panel of figure [3] It is symmetric around R = L/2 and thus invariant 
under the replacement R ^ L ~ R. The previous construction also entails that V{R) has to be extended 
periodically outside the range < R < L, although we have not performed that calculation explicitly. 
The level crossing at i? = L/2 is clearly visible, when the two Wilson states from figure [5] exchange their 
roles. The effective string tension is 



Cl+l 



_ dV{R) 



dR 



= o9 



(57) 



fl/L-s-O 



32 



33| or explicit quark wave functions 



and this agrees with eq. psp and calculations using path integrals 
35| . In the infinite volume limit L — oo, the potential is exactly linear, since all Wilson states other 



than I R) have energy of order ^{L) and thus decouple from the spectrum. 



V. YANG-MILLS THEORY IN (3 + 1) DIMENSIONS 

Let us finally address the most complex example, viz. non-Abelian Yang-Mills theory in (3+1) dimensions. 
We employ again the first method based on the standard gauge fixed Hamiltonian eq. (|24p . 

(ly) =tr(R^|e-^(^-^°)|R) S)A^J[A^]^R[A^]e-^'^"-^"^'^R[A^]. (58) 

Here, ^ is the Faddeev-Popov determinant and ^^^[A^] = ( A^ | C/(R) | 0) is the wave function of the 
Wilson state | R ) in Coulomb gauge. 

To proceed, we could take a model for the true Yang-Mills ground state in Coulomb gauge, | ) and 
compute ^Pj^JA^] by application of the Wilson line. Below, we will instead apply an approximation 
where only the vev. of the Hamiltonian (rather than the full evolution operator) is required. In such 
cases, it is expedient to perform the same unitary transformation eq. (|15p which transfers the Wilson line 
in the Hamiltonian, this time however with the gauge fixed version eq. (|24p of the initial Hamiltonian. 
Retracing the steps that led from eq. ([5]) to eq. (fT7|) with the parallel transporter ?7(R) now depending 
on A^ only, we find that a transversal electric field is induced which shifts the momentum operator in 
eq. ([24|) . Moreover, the transformation of the Coulomb term eq. ([25| replaces the full charge {p + pdyn) 
with (p + pdyn + Pind), whcre p is the rotated external charge eq. (fT9|) . pdyn = , n"*" ] and pi^d = 

— [ A"*- , c*" ] . The effective gauge fixed Hamiltonian is therefore 

i/*^'^ - H^Ij + Vc + AH, (59) 
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where the first piece is just the usual gauge-fixed Hamiltonian in the absence of external charges, 

,2 



rr fix 



J d\/-^IVi{^)JUi{^) + ±^ J d3x[B°(x)]^ 



d^x,y) ^-VSy„(x) / ■ F'"'(x,y)p^y„(y) 



The two additional pieces comprise a non-Abelian Coulomb potential 



(60) 



(61) 



and a new contribution which contains the coupling of the external charges p to the dynamical charge of 
the gluon, as well as the terms induced by the electric field due to the Wilson line, 

AH^^ J d3x{^-in^(x)/e^(x) + e^(x)n^(x)+[e^(x)]'} + 



J ^ PdynJ' F{p + Pind) + {p + Pind) F pdyn + 



+Pind F Pind + And F p + p F Pi, 



While a full calculation of (W) cannot be done with a Hamiltonian as complicated as eq. (ISH)) . it naturally 
lends itself to approximations based on models for the (uncharged) vacuum wave functional, since the 
whole effect of the Wilson line is already incorporated in the dynamics, i.e. in the transformed Hamiltonian 
H. In the following, we will sketch such an approximation based on a Gaussian model for the vacuum 



wave functional, as obtained recently by variational calculations 13l-l20|. 
In a first step, we use Jensen's inequality to obtain a lower bound for the Wilson loop (and thus an upper 
bound for the Wilson potential). 



'W^) ^ 



tr e 



-T{H''^-Eq) 



) > tr exp 







T{H^''-Eo) 







The advantage of the rhs is that only the vev of the Hamiltonian is required. In particular, the eigenvalue 

rfix 

'ym 



equation |0) = i?o | 0) entails that the YM Hamiltonian drops out and we are left with 



[W) > tT exp 







T {Vc + AH) 







(62) 



To compute the exponent on the rhs of this equation, we write out AH explicitly, 

^2 



0\Vc + AH\0] 



[0\Vc\0 



He" 



ah) 



(The explicit form of the remainder AH will be given below.) Next we neglect the commutators 



^]«0, 



[C/^(R), T°] « 



(63) 



where 1/^(11) is the Wilson line — > R built from the transversal vector field A-'-. The first equation 
means that the gauge field dependence of the Faddeev-Popov determinant is neglected, ^ = const, while 
the second equation discards the matrix structure of the non-Abelian parallel transporter, ^ 1. (This 
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is stronger than just neglecting path ordering). As a consequence, the external charge and the induced 
electric field also become independent of the gauge field, 

p°(x) ^ ^T° [(5(x) - 5(x - R)] 
[e^] 1 (x) ^ iT'^ I dvk Uk (x) (5(x - y) . (64) 
The relevant matrix element in the exponent of eq. ()62|) now reads 



(0|Vc + Aif|0)«^ j pFp+^^ J d'x[eii^)]' + (0\AH\0) (65) 

where we used global color and translation invariance to infer the structure 

(0|n,i(x)|0) =0, (0|F'^''(x,y)|0) =(5"''F(x-y). (66) 

The first term on the rhs of eq. (|65p is the non-Abelian Coulomb interaction, i.e. the remaining terms 
describe, within our approximation scheme, the difference between the (non-Abelian) Coulomb potential 
and the Wilson potential. We will now investigate these corrections: 



The third term on the rhs of eq. (p5)) involves the coupling of the various charges via the Coulomb kernel, 

Pdyn {p + Pind) F d + Pind F p + p F pind • (67) 



AH='- 
2 

Its vev is still too complex for a direct evaluation. To further simplify it, we make yet another approx- 
imation and assume that the Coulomb kernel, F, inside charge interactions can be replaced by its vev 
eq. (|66ll . In an obvious notation, this means 

(0|g''i^"''Q"'|0) w (0|Q"Q"'|0) • (0|F"''|0) (0|Q"Q'"|0) -^^ (68) 



for any two operators Q,Q' . We can now use the previous approximation eq. (j63p and consequently 
eq. iMl) to put ( I AH \ 0) in the form 

r,2 f ^ 

(Pdyn) Fp + (pdyn Pind) F + pF (pdyn) + (Pind Pdyn) F + (pind Pind) F + (pind) Fp + pF (pi„d) 

Since Coulomb gauge does not single out a color or space direction, global color and translation invariance 
tells us that 

(0|A^(x)|0)=0, (0|Af(x)n,^(x)|0)^5., 
so that (pdyn) — (Pind) = 0. Furthermore, we have the charge correlator 

( Pdyn Pind ) = ( Pind Pdyn ) = 6^ "^(x) ( | f"' A^ ^(x) A,^ ^(x) H,^ "^(x) | ) . 

(For the first equality, we have to commute 11 to the right and use (A) = 0.) The vev is a single 
component of a color and space vector, and again global color and (spatial) rotation invariance tells us 
that it has to vanish. We are thus left with only a single term 

( I Ai? I ) = ^ r"'^ j (x, y) ei'i^) (x) F(x - y) • ( | ^(x) ^ (y) | ) . 
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Using global color and translation invariance one last time, the gluon propagator must take the form 

(0 I 4^^(x) Ai^y) I 0) = i.fc(x) I?(|x - y|) (69) 
with an unknown scalar function D{r). Thus eventually, with /"'"^ f"''^'^ — Nc S^'^, 

2 

9 



:0\AH\0] 



e^(x)F(x-y)i?(x-y)e^(y). 



(70) 



The first term in eq. (1651) is the non-Abelian Coulomb interaction. It can easily be computed within our 
approximation scheme, since the simplified form eq. (|64p for the external charge implies 



We observe that this expression can be combined with the second term in eq. ([65 



;0|^c|0) = |- / pFp 



■Co 



F{R) - F{Q) 



(71) 



with the quadratic Casimir C2{Nc) — -~J2a T°-T°-. Putting everything together, we arrive at 

(PF)«exp| -T^iVc / d^x,y)ei{^) [iV^i 5(x - y) + F(x - y) Z?(x - y)] e,i(y) I x 



X exp 



F{R) - F(0) 



(72) 



To estimate the contribution from the first line, we require some dynamical input on the gluon propagator 
D and the Coulomb kernel F . (Within our approximation, this is the only place where the actual form of 
the vacuum wave functional enters.) The exact analytical expressions are unknown, but we have reliable 

which are also in decent agreement with high precision 



data from recent variational calculations 
lattice measurements . The latter give 

1 



D{k)^ 



F{k) ~ Stt 



2w(fc) ' 



fc4 ' 



w(fc) « v'PTmvP 



fc ^ 0, 



(73) 



where M k, 880 MeV is the Gribov mass and ac is the Coulomb string tension, which is about 2 ... 3 
times the Wilson string tension. 

With this input, it is shown in appendix |B] that the integral in the first line of eq. ([7^ yields a Coulomb 
type of potential which decays as 1 /R at large distances and is thus subleading as compared to the second 
line (which gives a linearly rising potential). Thus, 



const • exp | - T C2 5^ F(i?) + subleading | 
and the Wilson string tension equals the Coulomb string tension. 



(74) 



This simple result should, however, be taken with some caution: In deriving it, we had to make a number 
of approximations, 

1. Jensen's inequality eq. (p^ 

2. the commutator approximations eq. (|63p 
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3. the factorization eq. ([68ll 

4. the variational solutions for the Green functions D{r) and -F(x) 

The corresponding errors are not all under good control: While lattice calculations indicate that the 
variational solutions in item ^4 are close to the exact results, Jensen's inequality entails that our 
potential is, a priori, only an upper bound to the true Wilson potential, and it is hard to predict the 
quality of that bound. Even more severe is the factorization #3, since it causes most of the couplings 
between the various charges to vanish. In particular, we have lost the coupling of the external charges 
at the end points of the Wilson lines with the dynamical charges of the gluon, as described by the 
contribution 



It is tempting to expect that eq. (I75p will screen the Coulomb interaction when treated non-perturbatively, 
and that this effect provides the major reduction of the Coulomb string tension to the Wilsonian one. 
A detailed study of the screening of external color charges by dynamical gluons is subject to ongoing 
research. 



In this paper, we have investigated the temporal Wilson loop in the Hamiltonian formulation of Yang- 
Mills theory. Using a unitary transformation we have demonstrated that the effect of the Wilson lines 
is to introduce a new non-Abelian (field-dependent) electric field in the Hamiltonian. This formulation 
leads to the correct result in some simple cases where the exact gauge-invariant ground state is known. 

Without the exact ground state at hand, we need to fix the gauge, for which Coulomb gauge is a 
particularly useful choice: Since the only constraint on physical states is Gauss' law which can be resolved 
exactly, any normalizable Ansatz for the ground state wave functional is admissible and the formulation is 
ideal for variational approaches. We have shown how to treat the exactly solvable cases in Coulomb gauge, 
and applied these methods to the physically important case of Yang-Mills theory in (3 + 1) dimensions 
using the optimal Gaussian wave functional from recent variational approaches as input. Unfortunately, 
the dynamics of the system was still too complex so that additional approximations were necessary to 
find an analytical result. Within these approximation, the Wilson string tension was found to match 
the Coulomb string tension, cr « CTc. Finally, we have discussed the qualitative effect of the neglected 
couplings in the Hamiltonian, which are expected to suppress the Wilson string tension towards the value 
(t/(Tc — 0.3 ... 0.5 favored by lattice simulations. A closer investigation of those missing couplings is 
currently underway. 




(75) 



In perturbation theory, one can show explicitly that this term reduces the Coulomb potential j pF P- 



VI. CONCLUSIONS 
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Appendix A: The Abelian Coulomb potential in spherical coordinates 



It is convenient to use spherical coordinates x = (r, a;), R = (i?, R) and express the (5- function by 

<5(x-.R) =^±-^J2YLii)yimiR) ■ (Al) 

l,m 

With this representation we find with ((20)) 

1 1 oo 



X 







^ E E / (^)^'™' (^)^" (^)^'™' (^) • ( A2) 

Performing the integrals over r, s' and Qx by means of the orthonormahty of Yi^n{x) we obtain 



2 d^xe(x)^ = ^- / ^5i:^^(0), (A3) 



(2) 

where is defined by 



Ira 



(2') 

Note that 5^ (il) is dimensionless. From the final integral 



^ 1 



1 

= -l + cx) (A5) 

s=0 

we ignore the divergent contribution from the lower integration limit, which represents the self-energy 

(2'] 

contribution of the two charges. The divergent factor (5<t (0) has a different origin: it is due to the 
infinitely thin Wilson line between the charges and gives rise to the renormalization of the Wilson loop 
operator, cf. eq. ([33]). In the present context, this amounts to the renormalization of the bare charge. 



9'5'<\^) = ^^ (A6) 



so that the Coulomb potential from eq. (jASP becomes eventually 

V{R)^'^js.e{^f = -^. (A7) 

Our result (|A3p disagrees with the similar calculation in ref. |25|], where the divergent factor 

is replaced by 1/ (An) , so that no renormalization seems necessary (except for the subtraction of the 

self-energy) . 

To relate the present calculation to our reasoning in section III Al of the main text, it should be noted 
that eq. ()Aip implies s = r/R = |x|/i? and the lower bound s = in eq. (jA5| should therefore be 
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Smin — |x|min/fi, witli |x|inin & minimal distance of the two charges. If we put a cutoff Smin = a/i? for 
the lower integration limit in eq. (|A5|) . the (unrenormalized) potential becomes 

Clearly, the self-energy of the charges corresponds to the second term in the parenthesis. However, even 
with a minimal distance cutoff a > 0, there is still a divergent prefactor coming from the infinitely thin 
Wilson line. To see this, recall that e(x) has support on the Wilson line only, so that x = R if the Wilson 
line is infinitely thin, and hence the argument of the angular delta-function in eq. (jA3p vanishes. The 
divergent prefactor is thus well understood as a typical OPE divergence for the thin Wilson line which 
involves field operators A(x) at arbitrarily close positions. It must be absorbed in the multiplicative 
renormalization Zw of the Wilson loop operator, cf. the discussion following eq. (|33[) in the main text. 



Appendix B: Contribution of to the Coulomb potential 

To estimate the role of we calculate the vacuum expectation value of the induced Coulomb term 

■,2 



Using the approximation 



Ai7(i) = |. / A • e^Fk ■ . (Bl) 



{AFA)^(^AA)(F)^^ (B2) 



and ignoring the path ordering in e^(x), i.e. Tjyx) ~ T^, we obtain 

2 ^ K 

(Ai/^i)) ^C2^Jdx,J dyjUj{x)u:~'{x,y)F{x,y) , (B3) 



where C2 1 = —Ta Ta is the quadratic Casimir operator in the color representation of the Wilson loop. 
Putting the path along the z-axis from to i? 

Xi = disRs , j/i = SaRt (B4) 

we find after Fourier transformation and using spherical coordinates in momentum space 

^00 00 11 



1 1 

X / dz (1 - z2)g»2feii(.-t) f dz' e^'^'P^'-'-'K (B5) 



-1 -1 

We are interested in the large R behavior of this quantity. Assuming the usual IR-behavior of the Green 
functions, 

uj{k) ~ 1 , F{k) ^ 1 (B6) 
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as fc — ;> 0, and furthermore introducing dimensionless variables 

k = kR, p = pR, (B7) 

on finds 

{AhP) - 4 as i? ^ oo . (B8) 

-ft 

This quantity is therefore subleading at large distances and hence irrelevant for the screening of the 
Coulomb string tension. 
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